function img_out = he(img_in)
img_in = im2double(img_in);
r = histeq(img_in(:,:,1));
g = histeq(img_in(:,:,2));
b = histeq(img_in(:,:,3));
img_out = im2uint8(cat(3,r,g,b));
end